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1. Introduction 

The twisted mass formalism ^ has proven itself to be a powerful tool to extract precision 
physics in the pseudo-scalar |j3|, 0] and baiyon sector [Q]. In this work we look at one of the 
simplest resonances the p meson. The additional complication for the p meson is that it decays via 
the strong interaction to two pions. In this project we aim to study the effect of the decay of the p 
mass and decay constant via the chiral extrapolation formulae. As an extension of the study of the 
decay constants of the light vector meson we also studied the decay constants of the and r\c 
mesons. 

The twisted Wilson action with the tree level Symanzik gauge action was used (see [Q] for 
details). In this analysis we use the rif = 2 unquenched data sets at j8=3.9 24^ 48 and j8 = 4.05 
32^ 64 with lattice spacings determined from to be 0.0855 fm and 0.0667 fm respectively [^]. 
These lattice spacings were consistent with those from the mass of the nucleon |^. 



2. Mass of the p meson 

The masses and decay constants for the light vector meson were obtained by fitting an order 4 
smearing matrix to a factorising fit model [Q. The correlators were generated with the "one-end" 
trick to improve the signal to noise ratio. We only consider unitary data. Vladikas et al. studies 
an independent set of vector meson correlators with partial quenching on the same set of gauge 
configurations. 

In this section we discuss the chiral extrapolation of the mass of the p meson. Bruns and 
MeiBner have published a "new" chiral extrapolation formulae for the p meson, based on a 
modified MS regulator. 

Mp=M''p+c,Ml + C2Ml + C3Ml\n{j^J (2.1) 

It is important to check that size of c, from the fits to the lattice data is consistent with other 
estimates from phenomenology. Bruns and MeiBner [j8|] claim that phenomenology prefers | c,- | < 

3. Leinweber et al. |^] claim to know sign of C2 (negative) and magnitude from (opn coupling. 



The original effective field theory calculation Jenkins, et al. [ |10| ] had C3=0. There are too many 
parameters to determine from our data, so we use an augmented with the above constraints built 



in [11] 



Xaug X ~l~ 2^ 0.2 

We also tried the fit model suggested by Leinweber et al. [^]. The formulae doesn't obey the 
power counting of effective field theory, but does include an explicit term for the decay of the p to 
two mesons. Also some coefficients are fixed from experiment in the fit model. 

Although we have data at two lattice spacings we only fit the j3 = 3.9 24^ 48 data sets. Un- 
fortunately the p correlators at j8 = 4.05 and the 32^ 64 jS = 3.9 data were too noisy to be useful. 
See [ p^ for a comparison between our results for the mass of the p meson and other lattice collab- 
orations. We are investigating various measurements techniques, such as colour diluted stochastic 
sources, to reduce the statistical errors. 
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Figure 1: Chiral extrapolation of the mass of the p meson 

The lattice data are shown in figure [T[ A fit linear in the square of the pion mass gives nip 
= 867(29) MeV (j8 =3.9 data), compared to the experimental value of 770 MeV. Naively the p 
data prefers a different lattice spacing to that from fji, but we think that the error is due to missing 
chiral corrections. We fit the Bruns-MeiBner model to the data with constraints on the C2 and C3 
coefficients and obtain the preliminary result is nip = 807(8) MeV. For the Adelaide fit method 
for the p meson, we obtain mp = 847(72) MeV. The statistical errors need to be reduced for a 
detailed comparison between the lattice data and the results of effective field theory. 

3. Decay constants of the p meson 



Light cone sum rules require the transverse decay constant of the p meson [|13|, [14|] for the 



extraction of |^ from the B — > py and B — > K*y decays, and other semileptonic decays of the B 



meson [15] 



In the continuum the decay constant of the p meson is defined via 

(0 [ \if{x)Y^Y{x) \ p) = mpfpE^ (3.1) 
The transverse decay constant {fy{}l)) of the p meson is defined by 

(0 I VAa^vV^ I P) = ifv {l^){p^i£v-Pv£^i) (3.2) 
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i8 


Za 


Zr(At = l/a) 


3.90 
4.05 


0.771(4) 
0.785(6) 


0.769(4) 
0.787(7) 



Table 1: Non-perturbative renormalisation factors used in this analysis 



Group 


Method 


fjil GeV) MeV 


/pMeV 


it 
fp 


Becirevic et al. | 
Braun et al. [ 20 
QCDSF 2005 [| 
RBC-UKQCD 


^5] 
1 

I|] 

,2] 


quenched 
quenched 
unquenched 
unquenched 


150(5) 
154(5) 
168(3) 
143(6) 


256(9) 


0.72(2)+2 
0.74(1) 
0.66(3) 
0.69(3) 


This work (preliminary) 


unquenched 


177(26) 


229(35) 


0.76(14) 



Table 2: Summary of results for transverse and leptonic decay constants of the p meson 



where a^v = '72[7^,7v]. The fy[pL) decay constant can not be determined from experiment. For 
the charmonium part of this analysis we studied the decay constant of the pseudoscalar meson. 

To renormalise the currents we use the results from the Rome-Southampton method, slightly 
updated from those reported by Dimopoulos et al. [ [l^ ] last year. We remind the reader that in 
twisted mass QCD the charged vector current renormalises with Za- The renormalisation factor 
of the tensor current depends on the scale, so we run to 2 GeV using the method described by 



Becirevic et al. [|15p. We wiU include the results from Gracey's three loop calculation at a later 
time 0. 

The chiral perturbation theory calculations for the chiral extrapolations of the leptonic vec- 
tor decay constants have been calculated [jlsl]. The corrections due loops start at ni^logm^ and 
mj^ [p^]. Unfortunately our light decay constant data is too noisy to look for chiral corrections, 
so we use simple linear fits in the quark mass. The chiral perturbation theory for tensor sources has 
been developed, but no loop calculations are available ||l9|]. Although by taking ratios of correlators 
the ratio of j- can be extracted directly []15[], we prefer to separately compute fp and fp, because 



fp is known from experiment (207 MeV) is thus a good validation test that the p to In decay is 
under control. In table |2| we present our preliminary results for the light vector decay constants, 
and compare to previous results in the literature. 



4. Decay constants in charmonium 



One way to understand charmonium production and decay is to use the NRQCD formalism, 
where non-perturbative information is encoded in a few matrix elements. NRQCD mostly produces 
a good description of experimental data. See [23] for a review and a discussion of the convergence 
of the velocity expansion in the charm region. However the results from Belle [24] for double 
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charmonium production 

^7/i^+T]J^ = 25.6 ±2.8 ±3.4 /Z^ (4.1) 
where =^ < 1, are much larger than the prediction from leading order NRQCD by Bodwin et 



al. [25] 



o[e+e- ^ J/w + ^c\nrqcd-^o) = 3-78 ± 1.26 fb (4.2) 

BaBar has a similar result to Belle's for this process [^]. Although further NRQCD calculations 
that include relativistic corrections and higher order perturbative corrections have increased the 



leading order NRQCD results, the final result is still lower than the experimental result [ 27 ] 



It is claimed that calculations that use light cone wave functions of charm mesons agree with 



the Belle result []28[]. However, Bodwin [ ]29| ] remarks that it is not clear that the model light cone 
wave functions used are close to true quarkonium wave functions. Lattice QCD calculations should 
be able to constrain some of the parameters of the light cone wave function of the charmonium 
mesons. In particular the decay constants are parameters of the light cone wave functions. We refer 
the reader to figure 6 in | |30| ] to see the effect of the decay constants on the production cross-section. 



See 1 3 1 ] for a recent review of developments in this field. 



As stressed by the FNAL and HPQCD collaborations []32[], the computation of the decay con- 
stant of the meson, that is known from experiment, is an important validation test for calcula- 
tions that compute decay constants that contain heavy quarks. 

The lattice calculation is essentially the same as for the light quarks in section ^. As normal 
the charmonium correlators are much less noisy than for the equivalent light quark correlators. 
We used charged pseudoscalar interpolating operators, and we do not include any disconnected 
diagrams in the analysis. We computed the leptonic and transverse decay constant of the 
meson and the decay constant of the T]c meson for the j3 values 3.9 and 4.05. We used three heavy 
quark masses that interpolated the mass of the charm quark. We tuned the charm mass by looking 
at the 7/ 1//^ meson. Only local-local correlators were used. 

The results for the three decay constants are plotted in figure ^. We quote preliminary results: 
fjl^ = 413(40) MeV, fj/^-tram =396(35) MeV, and /^^. =379(29) MeV, for the decay constants in 
the continuum limit. The errors on the final results are inflated because we only had data at two 
lattice spacings. 

The experimental value of fj/^= 41 1(7) MeV (see [|^] for example). There is a "sort of exper- 
iment result" frj^ = 335(75) MeV from B — > rjcK with factorisation assumption from CLEO [p4|]. 
We now compare to other lattice QCD calculations. The Jlab lattice group obtained [ p3| ] /y/^, = 
399(4) MeV and /r,^ = 429(4)(28) MeV from quenched QCD, at a single lattice spacing. Chiu and 
Hsieh obtained /,,^ = 438(5)(6) MeV from a quenched QCD calculation using the overlap fermion 



action [35], at a single lattice spacing. 



5. Conclusions 

We have described our lattice QCD calculations to compute the mass and decay constants 
of the rho meson. We are trying to carefully study the chiral extrapolation. Similar to the chiral 
pertubation theory studies with light pseudoscalar mesons, the errors on the masses and decay 
constants have to be very precise to see the effects of loop diagrams. 
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* J/\|/ trans 
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Figure 2: Decay constants as a function of square of lattice spacing 

We reviewed why all three decay constants of the charmonium system are important for the 
phenomenology of double charmonium production and we presented preliminary results for these 
decay constants. We note that the moments of the light cone wave functions could also be deter- 
mined in a similar manner to the calculation in the light quark sector p^. 
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